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ABSTRACT 

The transient atoll source 4U 1608-52 had been extensively observed by the Rossi 
X-ray Timing Explorer {RXTE) during its 1998 outburst. We analyse its X-ray spec- 
tra as a function of inferred accretion rate from both the 1998 outburst and from the 
non-outburst 1996 and 1998 data. Wc can fit all the spectra by a model in which 
seed photons from the neutron star surface are Comptonized in a boundary layer. The 
Comptonized emission illuminates the accretion disc surface, producing an ionized, 
relativistically broadened reflection signature, while the direct emission from the ac- 
cretion disc can also be seen. The evolution of the source can be explained if the main 
parameter driving the spectral evolution is the average mass accretion rate, which de- 
termines the truncation radius of the inner accretion disc. At low mass accretion rates, 
in the island state, the disc truncates before reaching the neutron star surface and the 
inner accretion flow/boundary layer is mostly optically thin. The disc emission is at 
too low a temperature to be observed in the RXTE spectra, but some of the seed 
photons from the neutron star can be seen directly through the mostly optically thin 
boundary layer. At higher mass accretion rates, in the banana state, the disc moves in 
and the boundary layer becomes much more optically thick so its temperature drops. 
The disc can then be seen directly, but the seed photons from the neutron star surface 
cannot as they are buried beneath the increasingly optically thick boundary layer. 
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1 INTRODUCTION 

Low Mass X-ray Binaries (LMXBs) with a neutron star pri- 
mary can be observationally divided into two main cate- 
gories, namely atolls and Z sources (Hasinger & van der Klis 
1989). This classification is based on changes in both spec- 
tral and timing properties as the source varies, and probably 
reflects difi'erences in both mass accretion rate, M, and mag- 
netic field, B. The Z sources have high luminosity (typically 
more than 50 per cent of the Eddington limit) and magnetic 
field {B > 10^ G) while the atolls have lower luminosity 
(generally less than 10 per cent of Eddington) and low mag- 
netic field {B ~ 10** G) (Hasinger & van der Klis 1989). 

Both atolls and Z-sources show spectral changes which 
form a well defined track in a colour-colour diagram. The 
atolls show hard, power-law spectra at low luminosities (is- 
land state), which are startlingly like the low/hard spec- 
tra of Galactic black hole binaries (e.g. Barret & Vedrenne 
1993; van Paradijs & van der Klis 1994; Barret et al. 2000). 
At higher luminosities the spectra are typically much softer 
(banana state; e.g. di Salvo et al. 2000; Piraino, Santangelo 
& Kaaret 2000). The Z sources always show soft spectra. 



but there are subtle spectral changes which show up as a 
Z-shaped track on a colour-colour diagram. Both types of 
sources move along their tracks on the colour-colour dia- 
gram on time-scale of hours to days, except the island state 
of atolls, where motion can take days or weeks. They do not 
jump between the track branches. Most of the X-ray spec- 
tral and timing parameters depend only on the position of 
a source in this diagram. This is usually parameterized by 
the curve length, 5, along the track. 

Until recently the atolls were thought to show a char- 
acteristic C- (or atoll) shaped track on the colour-colour 
diagram. However, mapping the full range of colour-colour 
changes for these LMXBs is difficult because the data avail- 
able from most of the sources show little variability. There 
are only a few transient LMXBs with enough observations 
covering wide range of luminosities, to follow the whole track 
traced out in the colour-colour diagram. In these transient 
systems the C-shaped path is only a subset of a larger track. 
At low luminosities, evolution within the island state forms 
an upper branch, turning the C into a Z (Gierlinski & Done 
2002; Muno, RemUIard & Chakrabarty 2002). 

While the colour-colour diagram is clearly useful as a 
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way of parameterizing the spectral changes, true spectral fit- 
ting gives much more information on the underlying phys- 
ical changes in the source emission. Here we look at the 
detailed spectral shape of the transient atoll 4U 1608-52 as 
a function of its position on the colour-colour diagram, and 
show for the first time that all the spectra can be fit with a 
model in which a Comptonized boundary layer around the 
neutron star illuminates the accretion disc, giving rise to a 
relativistically smeared (and often highly ionized) reflected 
component and iron line. We show that the spectral changes 
mapped by the colour-colour diagram are consistent with a 
picture in which increasing average mass accretion rate in- 
creases the optical depth of the boundary layer and decreases 
the inner radius of the accretion disc. 



2 THE SOURCE 

4U 1608-52 is one of the transient neutron star LMXBs, with 
outbursts triggered by the disc instability which reoccur on 
time scales of ~ 80 days to 2 years (Lochner & Roussel- 
Dupre 1994). However, little is known about the system as 
the optical counterpart is very faint (Grindlay & Liller 1978; 
Wachter 1997). The distance is likewise poorly known. The 
current best estimate of 3.6 kpc is from observations of flux 
saturated (radius expansion, Eddington limited) type I X- 
ray bursts (e.g. Nakamura et al. 1989). The inclination is 
constrained only by the fact that there are no eclipses in the 
X-ray light curve. The neutron star spin may be indicated 
by the recently detected burst oscillation in RXTE at 620 
Hz (Muno at al. 2001) 

Hasinger and van der Klis (1989) classified 4U 1608- 
52 as an atoll source from its spectral and timing proper- 
ties. Mitsuda et al. (1989) observed with Tenma a continu- 
ous spectral transition from the banana high-luminosity, soft 
thermal spectrum to the island state with its low-intensity, 
hard power-law spectrum. They successfully fit these spec- 
tra by a two-component model consisting of the multicolour 
blackbody from the accretion disc and the Comptonized 
blackbody from the boundary layer. Higher energy obser- 
vations with BATSE show that the island state power law is 
cut-off above ~ 100 keV (Zhang et al. 1996). An iron fluores- 
cence line is detected in the Tenma X-ray spectra (Suzuki 
at al. 1984) with intensity weakly correlated with the source 
flux (Hirano et al. 1987). Ginga observations showed spec- 
tra in varying stages between the island and banana states, 
where the iron line (modelled including the Compton re- 
flected component) increased as the spectrum softened and 
the source moved towards the banana branch (Yoshida et 
al. 1993). 

The power density spectrum of 4U 1608-52 in the island 
state resembles that of the black holes in the hard state 
(Yoshida et al. 1993). It is characterized by a flat top noise 
(Py oc I/", with the slope s ~ 0) below and a turn-over 
above a certain frequency, Ub ~ 0.2-1 Hz. There is also a 
low-frequency quasi-periodic oscillation (QPO) around 2-9 
Hz. RXTE has seen a similar spectrum, where the high- 
frequency noise can be represented by the exponentially cut- 
off power law (Yu et al. 1997) . In the banana state the power 
spectrum consists of very low-frequency noise below ~ 1 Hz 
(power law with s ~ —1), and an exponentially cut-off power 
law at higher frequencies (Mendez et al. 1999). 



Twin kilohertz QPO's were detected by RXTE (van 
Paradijs et al. 1996; Berger at al. 1996; Mendez et al. 1998a) 
with varying peak separation (Mendez et al. 1998b). A third 
kilohertz QPO has also been detected (Jonker, Mendez & 
van der Klis 2000). The kHz QPOs are only found in a lim- 
ited range of the colour-colour diagram, near the transition 
between the island and banana states (Mendez et al. 1999). 



3 DATA ANALYSIS 

We have analyzed RXTE observations from the 1998 
February-April outburst (archival numbers P30062, P30188, 
and P30502) and from low- luminosity observations in 1996 
December (P10094) and 1998 June, August and September 
(P30419). All these data belong to the Proportional Counter 
Array (PCA) gain Epoch 3. A conservative prescription to 
characterise uncertainties in the instrument response is to 
add 2 per cent systematics to the data for bright sources (e.g. 
Cui et al. 1997), irrespective of which detector/layer combi- 
nation and gain epoch are used. However, inspection of data 
from the Crab shows that in Epoch 3, when using only top 
layer data from detectors and 1, smaller systematics of 0.5 
per cent are more appropriate (Wilson & Done 2001). We 
have repeated this Crab analysis, and found that even using 
a single absorbed (flxed at Nh = 3.2 x lO'^^ cm~^; Massaro 
et al. 2000) power law gives a reduced ~ 1-1 with 0.5 per 
cent systematic errors. This clearly shows that there are no 
residuals in the response (e.g. features from the Xenon edge 
at 4.7 keV, or from the rapid decrease in effective area at 
low energies) of the top layer of detectors and 1 at this 
level at this time. This is not the case for all layers and all 
detectors, where the reduced = 2.2 with 0.5 per cent sys- 
tematic errors and there is a clear feature in the residuals 
at low energies. Thus we use only data from Epoch 3, and 
from the top layer of detectors and 1, and apply 0.5 per 
cent systematic errors. However, we have checked that our 
conclusions are robust even to using 1 per cent systematics, 
although this gives unrealistically low values. 

First, we created PCA light-curves for each energy chan- 
nel in 128 s time bins. We excluded all the data affected 
by the type I bursts. This gave us 240 ks of PCA data. 
The light-curves were used to build a colour-colour diagram, 
deflning a soft colour as a ratio of 4-6.4 to 3-4 keV count 
rates, and a hard colour as a 9.7-16 over 6.4-9.7 keV ratio. 
These count rates are corrected for the slow variations in 
the PCA response (van Straaten et al. 2000; see Gierliriski 
& Done 2002 for technical details applied here). The dia- 
gram is shown in Fig. ^. 

Next, we extracted the PCA and High Energy X-ray 
Timing Experiment (HEXTE) spectra corresponding to dif- 
ferent positions in the colour-colour diagram. We have as- 
sumed that 4U 1608-52 shows a Z-shaped track on the di- 
agram, with inferred accretion rate increasing along the Z: 
from left to right on the upper branch, through the diagonal 
to the lower branch (Gierlinski & Done 2002; Muno et al. 
2002). The diagonal branch is sparsely covered in this data 
selection. We divided the diagram into eight regions (boxes), 
as shown in Fig. |l|, and averaged the spectra over each box. 
Neutron star LMXBs are known to vary in luminosity at 
given colours, creating the so-called parallel lines in colour- 
intensity diagrams (e.g. van der Klis 2001). Therefore, within 
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1.8 2 2.2 2.4 

Soft colour (4.0-6.4 keV)/(3, 0-4.0 keV) 

Figure 1. Colour-colour diagram of 4U 1608-52. Each point rep- 
resents 128 s of the data. Slanted crosses correspond to the obser- 
vations in 1996 December, when the source was at a rather low 
luminosity, while the squares correspond to the 1998 outburst 
data. Boxes show selection of the data for spectral analysis. Our 
choice of the boxes (and their order) follows a Z-shaped track in 
the colour-colour diagram, which is appropriate for the atolls, as 
it has been shown recently by Gierlinski & Done (2002). 



one box we average spectra of different intensities. Though 
not necessary identical, these spectra should be very similar, 
considering nearness of colours. In the classic terminology of 
the atoll sources boxes 1-3 belong to the island state, boxes 
5-8 to the banana, while the status of the box 4 cannot 
be established without joint spectral and timing analysis. 
To avoid accumulating spectral data from periods distant in 
time, we have created spectra 1 and 2 from 1996 data only 
(slanted crosses in Fig. l]), and spectra 3-8 from 1998 data 
only (filled squares in Fig. |l|). The PC A spectra are accom- 
panied by the simultaneous HEXTE spectra extracted from 
cluster 0. We will refer to these joint spectra as S1-S8. 

For spectral fitting we use the X-ray spectral fitting 
package xspec version 11 (Arnaud 1996). The error of each 
model parameter is given for a 90 per cent confidence in- 
terval. We use the 3-20 keV PCA data and 20-150 keV 
HEXTE data. The relative normalization of the PCA and 
HEXTE instruments is still uncertain, so we allow this to 
be an additional free parameter in all spectral fits. 



4 SPECTRAL MODEL 

Unlike black holes, the accreting neutron stars have both a 
surface and a magnetic field. Therefore, their X-ray spec- 
tra are expected to be more complex than these of the 
black holes. As well as emission from the accretion fiow 
(direct emission from an optically thick. X-ray illuminated 
disc and/or an optically thin inner fiow, corona or active re- 
gions; e.g. Czerny, Czerny & Grindlay 1986; Esin, McClin- 
tock & Narayan 1997; Beloborodov 1999) there can be emis- 



sion from the neutron star surface and the boundary layer 
(e.g. Popham & Sunyaev 2001). 

Observationally, the X-ray spectra of LMXBs in the soft 
state are well known to require at least two different spectral 
components, but the interpretation of these is not unique. 
The so-called Eastern model (Mitsuda et al. 1984, 1989) 
assumes that the soft component is from the disc, while the 
hard component is due to Comptonization in the boundary 
layer or the inner disc. In the Western model (White, Stella 
& Parmar 1988) the soft component is a single-temperature 
blackbody from the surface, while the hard component is 
a Comptonized emission from the disc. While it is likely 
that both disc and boundary layer are Comptonized at some 
level, it seems more probable that the boundary layer should 
have a higher temperature than the accretion disc (Popham 
& Sunyaev 2001). 

The hard state X-ray spectra of LMXBs are similar to 
those of the black hole candidates in their low/hard spectral 
state. They are dominated by a power law, with some con- 
tribution from an additional soft thermal component. The 
power law is usually interpreted as Comptonization of seed 
photons in hot, optically thin plasma (see e.g. Barret et al. 
2000). 

Thus, a physically motivated model which can poten- 
tially fit the hard component in all spectral states is thermal 
Comptonization. An optically thin, hot plasma would pro- 
duce a power- law spectrum (island state), while a plasma 
with intermediate or high optical depth and lower temper- 
ature can create a softer spectrum in which the high en- 
ergy rollover due to the electron temperature can be seen 
(banana state). We model the Comptonized spectrum us- 
ing an approximate solution of Kompaneets (1956) equation 
(Zdziarski, Johnson & Magdziarz 1996). The model (here- 
after THCOMP, not included in the standard distribution of 
xspec) is parameterized by the plasma electron tempera- 
ture, Te and optical depth Te (or equivalently the asymptotic 
spectral index), and blackbody seed photon temperature, 

^sccd • 

The seed photon energy is very important in the Comp- 
tonization model, since it gives a low-energy cut-off in the 
spectrum. For typical mass accretion rates onto a neutron 
star the seed photons from either the inner accretion disc 
or neutron star surface should have fcTgeed ~ 1 keV. This is 
plainly not far from the observed X-ray bandpass, so simple 
analytic approximations (e.g. the COMPST model or a cut- 
off power law) for the Compton scattered spectrum are not 
valid, since they extend towards lower energies as a uncut 
power law (see e.g. Done, Zycki & Smith 2002). The model 
THCOMP we use in this paper treats the seed photons low- 
energy cut-off properly. 

The boundary layer should illuminate the inner disc, 
producing a relativistically broadened iron line and associ- 
ated reflected continuum. While a broad iron line is often 
identified in atoll spectra, the reflected continuum which 
must accompany any line emission has only been fit to 
island state spectra. We calculate the Compton refiection 
of the THCOMP continuum spectrum using a model where 
both the continuum reflection and the line are calculated 
self-consistently for a given ionization state (Zycki, Done & 
Smith 1998). The reflected spectrum (continuum and line) 
is also relativistically smeared for a given inner disc radius, 
i?in. We approximate relativistic smearing by convolving the 
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Figure 2. PCA/HEXTE spectra of 4U 1608-52 along the atoll track, as selected from the data in Fig. |l|. The inferred accretion rate 
increases from SI to S8. The individual panels show the unfolded data and models, as well as the residuals. The HEXTE data have been 
rebinned for clarity. The model, consisting of the disc blackbody (the soft component), Comptonization (the hard component) and its 
reflection, is described in Section ^ The best-fitting parameters are shown in Fig. ^. 



reflected spectrum with a relativistic line profile (Fabian et 
al. 1989). 

First we demonstrate explicitly the need for a two com- 
ponent continuum in all the data sets. The thcomp contin- 
uum model and its reflection alone cannot give satisfactory 
fits for the spectra S1-S5, where > 120 at 80 or 79 de- 
grees of freedom. There is a strong soft excess in the resid- 
uals below ~ 5 keV, which is particularly marked in SI and 
S5 where > 200. We use a multicolour disc blackbody 
(diskbb in XSPEC; Mitsuda et al. 1984), to describe the soft 
component. This model is parameterized by the disc temper- 
ature, Tsoft, at the inner disc radius (Rin) and normalization 
A'soft oc -Rfjj. However, its shape is poorly constrained by the 
RXTE data so it can equally well be a single-temperature 
blackbody. The spectra S6-S8 give good fits without the soft 
component (x^ ~ 80, 88 and 75 at 78 d.o.f., for S6, S7 and 
S8, respectively), though even here inclusion of diskbb im- 
proves the fits significantly. The soft component in the case 



of the weakest improvement (for S6 where Ax^ = 11.6) is 
statistically significant at 99.7 per cent level. 

The Compton reflection (including the self-consistent 
iron line emission) is always significantly 
detected also. Using the two component continuum model 
(DiSKBB-f thcomp) without a refiected component gives an 
increase in x^ by more than 50 in all spectra (corresponding 
to an extremely high significance level). 

These results motivate us to choose a model consisting 
of thermal Comptonization, its Compton refiection and the 
disc blackbody as the model with which to fit all the spectra. 



5 RESULTS 

For spectral fitting we use the model described in the pre- 
vious section. We apply absorption corresponding to a fixed 
galactic hydrogen column of Nh = 1.5 x 10^^ cm~^ (Pen- 
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Table 1. Fitting results of the thermal Comptonization (with Compton reflection) plus disc blackbody model to the spectra S1-S8. The 

1/2 

bolometric, unabsorbed flux, -Fboli wa-s estimated from the model. Units are as follows: kT^^tt and feTgeed a.re in keV, N^^^^ is in km, ^ is 
in erg cm and -Fbol is in 10"^ erg cm~^ s~^. Brackets in fcTe column denote the fits in which To was fixed, and its lower limits were 
calculated additionally (no upper limits were found). Results from this table are visualized in Fig. H. A detailed description of the model 
used is in Section ^ The individual spectra are selected from the data as shown in Fig. |^. The spectra and the best-fitting models are 
shown in Fig. H. 
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ninx et al. 1989). In the island state spectra, S1-S3, the 
HEXTE statistics at higher energies are not good enough 
to constrain the electron temperature of the Comptonized 
component, T^. The high energy cutoflt begins to be seen 
in S4 and S5, although again the uncertainties on Te are 
large. All of these spectra provide only lower temperature 
limits. Therefore, we fit the spectra S1-S5 with kT^ fixed at 
50 keV, a value consistent with all these data sets, and then 
calculate the lower limits on the electron temperature. On 
the other hand, the spectra S6-S8 are able to constrain the 
electron temperature well. 

The data cannot constrain the inner disc radius from 
relativistic smearing of the reflected component. Spectral 
fitting gives only lower limits on _Rin, which are typically 
10-20i?g (where Rg = GM/c^) except for spectra S4 and S5 
which are consistent with Rin — 6Rg. We fix Rin = 20Rg in 
all the fits. 

We fit all the spectra with the above constraints, and 
show the spectra unfolded with the best-fitting models in 
Fig. 1^. The best-fitting model parameters are shown in Tab. 
land Fig. I 

We calculate the unabsorbed bolometric flux, Fboi, from 
the model. We cannot precisely estimate its uncertainties as 
this depends on the model used as well as on the statistics. 
However, fits with different soft component model (disc or 
blackbody) show consistently similar results, with the same 
trend: an increase from SI to S3, a drop at S4 and a further 
increase from S5 to SB. This supports the inferred direction 
of the increasing accretion rate on the colour-colour diagram, 
tracing out a Z-shaped track (see also Gierlinski & Done 
2002). 

With increasing inferred accretion rate the spectra be- 
come brighter, softer, and have a decreasing electron tem- 
perature as seen from the high energy cutoff. The change in 
the spectral shape between the island and banana states is 
clearly visible. Spectrum S4 bears more similarities to the 
banana state spectra than to island state, so it likely be- 
longs to the banana state. The curvature expected from the 
low energy cut-off in the Comptonized spectrum close to the 
seed photon energy is clearly seen in spectra S4-S8, where 
kTseed > 1 keV. The lower seed photon temperature de- 
rived for the island state means the low energy cutoff is less 
pronounced (not significantly detected) in spectra S1-S3. 



The electron temperature of the hard Comptonized 
component is consistent with a monotonic decline along 
the colour-colour track, while the optical depth is consis- 
tent with a corresponding monotonic increase. These two 
parameters are difficult to disentangle in data in which the 
high-energy rollover (at a few kTe) is not observed, as in the 
island state spectra. However, while the optical depth de- 
rived for a given temperature also depends on the assumed 
geometry, it is clear that the island state spectra have Te of 
order unity, while the banana branch has Te ^ 1. We derive 
similar optical depths and temperatures when we replace 
the approximate thcomp thermal Comptonization model 
with COMPPS Comptonization code (Poutanen & Svensson 
1996), which finds a numerical solution of the Comptoniza- 
tion problem explicitly considering successive scattering or- 
ders. 

Reflection (mostly driven by the detection of the self 
consistently produced iron line) is always significantly de- 
tected, although the inferred solid angle subtended by the re- 
flecting material is generally rather low, with (yi/2-n) ~ 0.1. 
The spectrum S4 has somewhat more reflection but this 
is not strongly statistically significant. Refiection is always 
strongly ionized, with ionization increasing in the higher lu- 
minosity banana spectra. The presence of reflection in all 
spectra (both in island and banana states) demonstrates 
that the broad, ionized iron line seen from these sources can 
indeed be produced by illumination of the accretion disc. 

The rather small reflected fraction in the island state 
is in conflict with Q,/2n « 0.5-1 found in 2-20 keV Ginga 
X-ray spectra of 4U 1608-52 by Yoshida et al. (1993) and 
Zdziarski, Lubiriski & Smith (1999). This is mostly due to 
the continuum model used. A power law and its reflected 
component does indeed give a larger Q/2tt, but these fits are 
statistically unacceptable, as well as physically inconsistent. 
Comptonization does not give rise to a power law at energies 
close to the seed photon energies, and the data also require a 
soft component as well as the Comptonized emission. Our re- 
sults are consistent (although with large uncertainties) with 
a reflection-spectral shape correlation for spectra S1-S4, as 
claimed by Zdziarski et al. (1999) although there is plainly 
no overall correlation that extends from the island state all 
the way through the banana branch. 

The observed soft component shows rather complex be- 
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Figure 3. Results from Tab. |l[ Units are as follows: kTsoCt and 

fcTseed are in keV, A'^Yft ''^ ? ^'^'^ Phol is 

in 10~^ erg cm~^ s~^. While fitting spectra SI— S5 Te was fixed 
at 50 keV, and its lower limits were calculated additionally (no 
upper limits were found). The symbols used here correspond to 
the symbols in Figs. ^ and ^ We note that the box number on the 
horizontal axis does not exactly correspond to the curve length 
in the colour-colour diagram, and in particular, there is a gap 
between boxes 3 and 4. 



haviour. Fig. q shows the soft component flux derived from 
the model given in Table 1, versus its temperature, Tgoft- 
There is a striking difference between the behaviour of the 
soft component in the island and banana states. The figure 
also shows the curve -Fsoft oc T^^f^, as expected if the disc 
geometry remained constant. The uncertainties are large, 
but while the banana state spectra could be consistent with 
this relation, the island state is clearly not. This can be also 
seen from the fitting results (Fig. where the normaliza- 
tion A'goft oc 7?f„ is fairly constant along the banana track 

The difference in the soft component between the is- 
land and banana can be further seen in Fig. ^. The seed 
photons temperature is systematically higher then the soft 
component temperature in the banana state. Though 90 per 
cent confidence errors for Tsocd and Tsott overlap (see Tab. 
|l]) , these two parameters are correlated and a fit with forced 
Tseod = Tsoft is significantly worse. In the island state the un- 
certainties on these temperatures are substantial but Tsccd 
is consistent with being equal to TsoSt, and the soft compo- 
nent can be the source of seed photons. This result does not 
change when we use a single-temperature blackbody instead 
of the multicolour disc. Conversely, in the banana state, the 
observed soft component is not the source of seed photons 
for the Comptonized component. 

In the island state, the observed seed photons (which 
are consistent with being the seed photons) have rather low 
normalization. If these were emitted from a disc then the 
inner disc radius can be estimated from the diskbb model 
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where D is the distance to the source, i is the inclination 
angle of the disc, /coi is the ratio of the colour to effective 
temperature (Shimura & Takahara 1995) and rj is the cor- 
rection factor for the inner torque-free boundary condition 
(Gierlihski et al. 1999: rj — 2.7 for Rin = 6Rg and less for 
higher Rin). With 77 = 2.7, D = 3.6 kpc, /coi = 1.8 and 
cosi = 0.5 we find Rin < 5 km for S1-S3. This is less than 
the expected neutron star radius of ~ 10 km, and suggests 
that though we fit the soft component by a disc model, it 
might not be a disc at all. Even if the colour temperature 
correction is as high as /coi ~ 2.7 (Merloni, Fabian & Ross 
2000), then S3 still yields _Rin < 5.5 km. There are of course 
more unknowns in Eq. (^, in particular the inclination an- 
gle, but Rin > 10 km requires i > 81°, which is excluded by 
the lack of X-ray eclipses. 

It is therefore unlikely that the soft component is the ac- 
cretion disc in the island state. Instead, we suggest it comes 
from the neutron star surface. A fit of a model including a 
single-temperature blackbody instead of a disc to the island 
state data yields an apparent radius of the blackbody of 4±2 
km. 

Conversely, in the banana states, the soft component is 
consistent with emission from the accretion disc, although 
these photons do not provide the seed photons for Compton 
scattering. The inner disc radius is consistent with remaining 
constant throughout the banana state. From the fit to the 
typical banana spectrum S6 we estimate Rm ~ 30 km, a 
value consistent with the disc truncated not very far from 
the neutron star. 

If the soft component is the neutron star surface in the 
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Figure 4. Soft component flux versus its temperature. The curve 
represents tlie best fit of a function F^^ft = -^"^sott ^° spectra 
S4— S8. The island state spectra S1-S3 are not consistent with this 
relation. 
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Figure 5. Seed photons temperature versus soft component tem- 
perature. The line shows were Tg^od = '^soft • In the banana state 
(S4-S8) 

'^seed is systematically higher than Tgof^ , therefore the 
soft component cannot be the source of seed photons. In the is- 
land state (SI— S3) the uncertainties are large, but it is possible 
that the soft component is the source of seed photons. 



spectra S1-S3 and the accretion disc in the spectra S4-S8, 
then the change of trend in Tsoft and A'^soft between S3 and 
S4 (see Fig. ^) is naturally explained, and does not contra- 
dict our assumption about the direction of increase in the 
accretion rate in the upper branch in the colour-colour dia- 
gram. We discuss a physical picture for the evolution of the 
soft and seed photons in Section 7. 



6 COMPARISON TO OTHER SOURCES 

The X-ray spectra of 4U 1608-52 analyzed here are entirely 
typical of atoll sources. Having a good coverage of a large 
luminosity variation, we can compare our results to previous 
observations of other atolls in various luminosity states. 

The island state is characterized by a hard, power-law 
spectrum, similar to the hard state of black hole binaries, 
though a bit softer. The presence of a soft component has 
been reported before in SAX J1808. 4-3658 (Gierliriski, Done 
& Barret 2002), 4U 1724-308 (Guainazzi et al. 1998), IE 
1724-3045 and SLX 1735-269 (Barret et al. 2000). In these 
sources, the radius of the soft component is not consis- 
tent with the accretion disc (except perhaps SLX 1735-269, 
though there is Galactic bulge emission component in the 
spectrum, which makes spectral modelling difficult; see Bar- 
ret et al. 2000). This is consistent with what we find in 4U 
1608-52 and suggestive that the soft component arises from 
the neutron star surface rather than from the accretion disc. 

The banana state X-ray spectra are soft, and here again 
the 4U 1608-52 is entirely similar to other atolls. A two- 
component model has been successfully applied to e.g. KS 
1731-260 (Narita, Grindlay & Barret 2001), GX 3-Hl, Ser 
X-1 (Oosterbroek et al. 2001) and MXB 1728-34 (di Salvo 
et al. 2000), the soft component being consistent with the 
disc. On the other hand, one must be very cautious here, 
since the soft component properties strongly depend on the 
Comptonization model applied, and inappropriate approxi- 
mations have often been used (see Done et al. 2001 for de- 
tails). 

The ambiguity in the spectral shape of the soft compo- 
nent has been reported before (e.g. Guainazzi et al. 1998; 
Barret et al. 2000). Interestingly, in MXB 1728-34, even 
BeppoSAX spectra which extend down to much lower en- 
ergies than the PCA cannot distinguish between a disc and 
blackbody spectrum for the soft component (see di Salvo et 
al. 2000). 

Compton reflection has been reported before in the is- 
land state, in e.g. SAX 1808.4-3658 (Gierlinski et al. 2002), 
4U 1728-34 (Narita et al. 2001), 4U 0614+091 (Piraino et 
al. 1999), GS 1826-238 and SLX 1735-269 (Barret et al. 
2000). However, this paper shows the first detection of re- 
flection in the banana state. Previous work has generally flt 
an ad hoc broad Fe Ka line, rather than including the self 
consistently produced reflected continuum e.g. in GX 3-1-1, 
Ser X-1 (Oosterbroek et al. 2001) and 4U 1728-34 (Piraino, 
Santangelo & Kaaret 2000; Di Salvo et al. 2000). 



7 DISCUSSION AND CONCLUSIONS 

We flt the spectra of the atoll 4U 1608-52 at all points along 
its track on the colour-colour diagram by a physically mo- 
tivated model, consisting of thermal Comptonization and 
its Compton reflection, together with a soft component. 
Both reflection (with the self consistently calculated iron 
line emission) and the soft component are always statisti- 
cally signiflcantly detected in the spectra. 

There is a dramatic change in spectral shape between 
the island and banana states, probably caused by an equally 
dramatic change in the accretion flow geometry. In the island 
state the observed soft component luminosity and temper- 
ature implies an emission region which is smaller than an 
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accretion disc around a neutron star. We suggest it arises 
from optically thick, thermal emission from surface of the 
neutron star. The hard component is consistent with Comp- 
tonization of these observed soft photons from the neutron 
star surface by hot electrons in an inner optically thin ac- 
cretion flow (or outer boundary layer). The optically thin 
Comptonizing medium has temperature of J; 20 keV. We do 
not see direct emission from the accretion disc in the PCA 
energy range, but we do see weak reflected emission imply- 
ing that the disc subtends a solid angle ~ 0.1 x 27r from the 
point of view of the hard X-ray source. 

In the banana branch spectra the Comptonized com- 
ponent is much softer, with temperature of ~ 5 keV and 
optical depth of ^ 1. Here the observed soft component is 
consistent with being from the disc, but the seed photons are 
not consistent with being the observed soft component, and 
again it is more likely that the seed photons are produced 
predominantly by the neutron star surface. This interpreta- 
tion favours the Eastern model of LMXBs on the banana 
branch. 

A scenario consistent with all of the above results is 
one in which the main parameter driving the spectral evo- 
lution is the average mass accretion rate, which determines 
the truncation radius of the inner accretion disc. At low 
mass accretion rates the optically thick disc does not ex- 
tend down to the neutron star surface, and the inner ac- 
cretion flow/boundary layer is mostly optically thin. The 
large radius of the optically thick disc means that its tem- 
perature and luminosity is low so it cannot be seen directly 
in the PCA spectra, and it subtends a rather small solid 
angle for reflection. Since the accretion flow is not very op- 
tically thick then some fraction ~ e~'^ ~ 0.3 of the emission 
from the neutron star surface can be seen directly, while 
the rest is Comptonized into the hard component. As the 
mass accretion rate increases, so the disc extends further in 
and the optical depth of the inner flow/boundary layer in- 
creases. Eventually the neutron star surface is shrouded by 
the increasingly optically thick flow, so that these seed pho- 
tons cannot be seen. The high optical depth also means that 
the boundary layer emission is then close to thermalizing, 
so its temperature drops, while the decreasing inner radius 
of the disc means that its emission starts to become more 
important in the PCA bandpass. 

In this scenario the island state/banana branch transi- 
tion is marked by the point at which the inner accretion flow 
collapses into a disc. The disc extends further in so that it 
is observed more easily as a soft component, softening the 
soft colour. But it also means that the boundary layer be- 
comes much more optically thick, so its temperature drops, 
softening the hard colour. 

It is plain that this can explain qualitatively the atoll 
source evolution. We will investigate the quantitative impli- 
cations of this in a later paper. 
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